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Abstract: We investigate the discrete j3 function of the 2-flavor SU(3) sextet model using 
the finite volume gradient flow scheme. Our results, using clover improved nHYP smeared 
Wilson fermions, follow the (non-universal) 4-loop MS perturbative predictions closely up to 
~ 5.5, the strongest coupling reached in our simulation. At strong couplings the results 
are in tension with a recently published work using the same gradient flow renormalization 
scheme with staggered fermions. Since these calculations define the discrete (3 function 
in the same continuum renormalization scheme, they should lead to the same continuum 
predictions, irrespective of the lattice fermion action. 

In order to test systematic effects in our computation we compare two different lat¬ 
tice operators, three different flow definitions, and two volume extrapolations. We find 
agreement among these different approaches in the continuum limit when the gradient flow 
parameter c > 0.35. Considering the potential phenomenological impact of this model, 
it is important to understand the origin of the disagreement between our work and the 
staggered fermion results. 
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1 Introduction 

The 2 flavor SU(3) sextet model has received considerable attention lately as candidate to 
describe the Standard Model Higgs boson as a composite particle [1-8] . Lattice simulations 
show that the model has a light O'*"*" scalar state that is well separated from the rest of 
the spectrum [9]. Quantum corrections could lower the scalar mass to match the Higgs 
boson, while the vector meson is predicted in the 2TeV range [10, 11]. If the model is 
chirally broken, it has 3 Goldstone bosons as needed for electroweak symmetry breaking. If 
it is conformal, an additional 4-fermion interaction might drive it into the chirally broken 
regime [12]. In either case the 2-flavor sextet model is a promising strongly coupled system 
that deserves detailed non-perturbative investigations. 

The infrared properties of the model are not yet settled. 2-loop perturbation theory 
predicts an infrared fixed point (IRFP) around ~ 10.58. The 3- and 4-loop (non- 
universal) MS predictions for the IRFP are lower, around ~ 6.28 and 5.90 respectively, 
while the mass anomalous dimension in all cases is 0(1) [13]. In a series of lattice cal¬ 
culations that used different improved Wilson fermions the authors of [14-16] found that 
the renormalization group (RG) /3 function of the Schrodinger functional gauge coupling 
became smaller than the 2-loop perturbative prediction for g^ > 4.0 and possibly developed 
a conformal IRFP around g^ ~ 6.0. Ref. [16] predicted the mass anomalous dimension of 
this fixed point around 'jm ^ 0.45. However these early works were not able to take the 
proper continuum limit. 
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There are several recent studies that use staggered fermions. In order to describe 2 
flavors the fermions are rooted. Rooting is not expected to be a problem when the proper 
continuum limit (/g —?• 0 is taken at a chirally broken system, though the procedure is 
less understood in conformal or walking theories near an IRFP. Finite temperature studies 
could distinguish conformal and chirally broken behavior but, as the results of refs. [17-19] 
demonstrate, the phase diagram of the 2-flavor sextet model is not easy to characterize. 
While the results of refs. [17, 18] appear to be consistent with chiral symmetry breaking, the 
latest results reported in ref. [19] favor the conformal scenario. The LatHC collaboration 
has been pursuing comprehensive, large-scale studies of this model for several years. Their 
results, presented in a series of papers, support chiral symmetry breaking [7-9, 11]. In 
ref. [20] the collaboration investigated the step scaling function, and found that for gauge 
couplings < 6.5 it is monotonic and shows no indication of an IRFP. 

In this paper we describe our investigation of the step scaling function using Wilson 
fermions. We use the same gradient flow renormalized coupling scheme as ref. [20] but find 
that the continuum limit extrapolated step scaling function significantly deviates from the 
2-loop perturbative prediction in the ~ 4.0 — 6.0 region, where it follows quite closely 
the (non-universal) 4-loop MS prediction. We consider various systematic effects that 
could influence our results. We compare gradient flow couplings defined with two different 
operators, the plaquette and clover forms. We also compare several different gradient flow 
schemes based on the t-shifted gradient flow coupling [21]. These gradient flow schemes 
predict the same continuum limit but have different lattice artifacts at finite bare coupling. 
We also vary the lattice volumes used in the continuum extrapolation. We find that if 
the gradient flow parameter c > 0.35, different variations predict the same continuum 
extrapolated step scaling function. Considering the phenomenological importance of this 
system, it is imperative to understand and resolve the more than 2>a tension between 
ref. [20] and our results. 

In the next section we briefly describe the calculation of the step scaling function using 
the gradient flow method. In section 3 we outline the lattice action used in this study, 
discuss the phase diagram, and problems related to the tuning of the fermion hopping 
parameter to the critical surface. Section 4 outlines the analysis that leads to our result 
while in section 5 we discuss cut-off effects. Section 6 concludes this study. 

2 The step scaling function with gradient flow coupling 

The gradient flow is a continuous and invertible transformation [22, 23] that has been 
used in a wide variety of applications. (For recent reviews see [24, 25].) The renormalized 
gradient flow coupling [26] is defined at energy scale /r as 

9'gf(/^) = 1 (2-1) 

where E{t) = —^ReTr [G{t)] is the energy density after flow time t and g, = l/v^- 
The normalization A/" is set by requiring that gQpif^) matches the continuum MS coupling 
at tree level. 
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At finite bare coupling the gradient flow renormalized coupling ^Qp has cut-off effects 
that depend on the lattice action, the gradient flow transformation, and the lattice operator 
used in defining the energy density E{t). While changing the lattice action is usually 
not practical, considering several flow transformations and/or lattice operators helps to 
investigate and reduce cut-off effects. In the present study we compare two lattice operators, 
the plaquette and the clover operators, to discretize the energy density [26]. We also 
consider several definitions of the gradient flow coupling based on the t-shift approach 
introduced in ref. [21]. The t-shifted gradient flow coupling is defined as 

~2 / ^ 2 / ^{E{t + Toa‘^)) 

S-GFlAi; “) = S-GFlAi; a) "-— (2.2) 

with |to| <C t/a^. The tq t-shift can be either positive or negative. In the continuum limit 
Too^ —)• 0 and gQp{fi) = i?gf(/^)- lowest order the effect of the t-shift can be described 
as an O(a^) correction to the coupling ^p(^;a); therefore, by optimizing tq the leading 
C>(a^) cut-off effects can be removed. In previous studies of the 4-, 8- and 12-flavor SU(3) 
systems [21, 27] the optimal tq depended only weakly on ^p(/r) and simply setting it 
to a constant value was sufficient to remove most observable lattice artifacts throughout 
the ranges of couplings explored in each case. Alternatively, by considering various values 
of the t-shift one can identify not only the leading O(a^) but the remaining higher order 
cut-off corrections as well [28]. Since the gradient flow is evaluated through numerical 
integration, replacing gQp t>y shifting t ^ t + tq can easily be incorporated in the 

analysis, without recalculating the flow. 

In step-scaling analyses one connects the energy scale to the lattice volume by fixing 
the ratio c = \^/L, as described in Refs. [29-31]. Different choices of c define different 
renormalization schemes, predicting different discrete /? functions in the continuum limit. 
If periodic boundary conditions (BCs) are used for the gauge fields, these /3 functions are 
only one-loop universal [29]. Ref. [32] proposed to use a modified renormalized coupling 


to perturbatively correct for cutoff effects. Here the function C{L,c) is a four-dimensional 
finite-volume sum in lattice perturbation theory, which depends on the gauge action, flow 
and operator. Its tree level expression is given in ref. [32]. Since we use periodic BCs for 
the gauge fields, the correction C{L,c) also includes a term that accounts for the zero-mode 
contributions [29]. 

We combine the t-shifted coupling of eq. (2.2) with of eq. (2.3) and use the resulting 
^ gradient flow running coupling to investigate the discrete /3 function corresponding to 
scale change s. 


Ps{9l\L) 


g^{sL;a) -g^{L;a) 
log(s2) 


(2.4) 


Our final results for the continuum discrete /3 function Ps{g^) = ^^^{a/L)^o Psid^, L) are 
then obtained by extrapolating (a/L) to 0. If c = y/^/L is fixed, the different values 
of To and different energy density operators should all predict the same Psig^) in the 
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continuum limit. Any deviation signals cut-off effects that are not removed by the (a/L) —)• 
0 extrapolation. These systematical errors can be controlled either by increasing the lattice 
volumes used in the continuum extrapolation or by increasing the parameter c [31]. 


3 Numerical details 
3.1 Lattice action 

Wilson fermions are often plagued by bulk first order phase transitions that prevent simu¬ 
lations to reach the strong coupling regime. Improving the action by smearing the gauge 
fields in the fermion determinant or by including additional gauge terms can help, as 
was illustrated in refs. [14-16], whose authors investigated the step scaling function using 
the Schrodinger functional running coupling of 2-flavor sextet fermions. Using thin link 
Wilson-clover fermion ref. [14] could probe the system in the range of < 2.4. Improving 
the action by using nHYP smeared fermions in ref. [15] extended the range to < 5.5. 
Including a smeared plaquette term in the gauge action pushed the first order phase tran¬ 
sition to even stronger coupling, where yet an other issue arises [16, 33]. On rough gauge 
conhgurations zero modes in the U(3) projection of the fermion smearing transformation 
can occur, a known problem in projected smearing transformations like nHYP or HISQ. 
To solve this problem ref. [34] added a pure gauge term suppressing the zero modes of 
the nHYP smearing to the action. This dislocation suppressing term, S'ndSj significantly 
reduces the occurrence of zero modes and allows simulations deep in the strong coupling 
regime. 

We have tested several versions of the nHYP smeared clover-Wilson fermion action 
for this project. We considered reduced nHYP fermions where the smearing coefficients 
(oi = 0.5 ,02 = 0.5 ,03 = 0.4) are chosen instead of the optimal original parameters 
(oi = 0.75,02 = 0 . 6,03 = 0-3). While the reduced coefficients do not provide the same 
suppression of the ultraviolet fluctuations as the original ones, they give significant im¬ 
provement [35]. At the same time reduced nHYP smearing cannot have zero modes in the 
smearing projection, making numerical simulations stable even without the S'nds term. We 
tried reduced nHYP smeared fermions with a simple plaquette gauge action, with a gauge 
action that included a smeared plaquette [16], and also tested a plaquette gauge action 
with a negative adjoint plaquette term [35]. While numerical simulations were stable, all 
three actions were limited by the occurrence of first order phase transitions that prevented 
simulations to reach couplings g^ > 4. 

Adding an NDS term to the action proved to be significantly better than any other 
choices. In this work we use nHYP smeared fermions with original smearing coefficients 
(ai = 0.75,02 = 0 . 6,03 = 0.3), plaquette gauge action and an NDS term with coefficient 
7 = (71 = 72 = 73 =)0.075. Perturbatively the NDS term shifts the gauge coupling as 

1 / 37.1 , , , 

“2 ^ Yv“ “2 + “ 3 ) (3.1) 

g^ 2Nc Nc 3 
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where /? is the coefficient of the plaquette gauge term. With our choice of parameters this 
shift amounts to 


1 _ +0.1725 

7o~ 2W 


(3.2) 


a fairly small effect. We found that simulations with this action are stable even at f3 = 0.5, 
though the step size of the molecular dynamics evolution has to be decreased significantly 
at strong bare couplings. 


3.2 Tuning to the critical surface 

Since our main goal is to determine the step scaling function, we generated configurations 
on symmetric volumes with Lja = 8 — 28. We chose standard boundary conditions 
(periodic for the gauge fields and periodic in space, antiperiodic in time for the fermions) 
and tuned the fermion hopping parameter n to the critical surface. 

The critical Kcr(/3) line in the (/3, k) plane is defined through the vanishing of the quark 
mass niq. At zero temperature, the chiral limit of the Wilson type fermion action is taken 
along this line. We define niq through the axial Ward identity (AWI), which relates the 
pseudoscalar density /2)Tp to the axial current A“ = ' 07 ^ 75 (r“/ 2 )' 0 : 

dt t)0“) = 2m, (P“(x, t)0“) , (3.3) 

X X 

where 0“ is a pseudoscalar Gaussian source and the summation is over the three spatial 
coordinates. To obtain the pseudoscalar meson mass (mps = Eq) and the quark mass m,, 
we fit the pseudoscalar (PP) and the axial vector (AP) correlators to the functional form 


N-l 




n=0 

N-l 


^(A“(x,f)0“) = Y,Cn,PpCn 


,AP 


n=0 


g ^nt _ g En{T t) 


(3.4a) 

(3.4b) 


From eqs. (3.3) and (3.4), we obtain the quark mass 

Eq C'o,ap 
= -^r —• 

2 Fo,PP 


(3.5) 


Figure 1 shows the «„(/?) line and Table 1 lists the corresponding values, together 
with the quark masses m, on 24^ volumes. We found that tuning to the critical line has 
to be done with increasing precision at strong coupling and on large volumes. Conformal 
systems in finite volume and at zero fermion mass are volume squeezed, the infrared cut-off 
is provided by the volume, oc 1/L. Increasing the mass does not have a strong effect as 
long as the system stays volume squeezed, m, <C 1/T. However at larger m, the role of 
the volume diminishes, the system moves into a mass-deformed regime. The spectrum in 
the mass deformed regime is considerably different than in the volume squeezed one. This 
phenomena was clearly observed in the finite size scaling analysis of the = 12 SU(3) 
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Figure 1. Critical kappa line determined from niq = 0 in eq. (3.5). 
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K 

mq 

f3 

K 

mq 

1.5 

0.14150 

0.0012(3) 

4.5 

0.13020 

0.0020(2) 

2.0 

0.13987 

0.0007(3) 

5.0 

0.12910 

-0.0007(1) 

2.5 

0.13800 

0.0018(3) 

5.5 

0.12830 

-0.0008(1) 

3.0 

0.13603 

-0.0005(3) 

6.0 

0.12763 

0.0028(1) 

3.5 

4.0 

0.13380 

0.13190 

0.0040(3) 

0.0004(2) 

6.5 

0.12727 

0.0016(1) 


Table 1. Quark mass values on 24^ volumes at k hopping parameter values used in this work. 


system [36]. Chirally broken but volume squeezed systems could show similar behavior [37, 
38]. 

In our simulations we observed that incorrect tuning of the fermion mass in either 
direction can lead to a sudden change in the gradient flow coupling , leading to severe 
lattice artifacts. In hgure 2 we show the gradient flow coupling as the function of c = y/^/L 
at several values of k on IG"^, /? = 1.5 configurations, both below and above Her ~ 0.14150. 
The K values in figure 2 are equally spaced. The sudden increase between k = 0.14120 and 
0.14100 signals that the system moved from the volume squeezed to the mass deformed 
regime. The qualitative change between k = 0.14160 and 0.14180 indicates the same 
phenomena, but now in the negative mass region. 

The transition between the volume squeezed and mass deformed regimes depends on 
the volume. Figure 3 shows the gradient flow coupling on L/a = 14 — 28 volumes. The left 
panel corresponds to (3 = 1.5, k = 0.14150, the right to /? = 3.5, k = 0.13380. We chose 
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1 ' I ' r 

L = 16; (3 = 1.5; % ^ = 0.15 



Figure 2. The gradient flow coupling as the function of c = v^/T on 16^ volume and /3 = 1.5 
at several, equally spaced, k values. Kcr ~ 0.14150 at this gauge coupling. 




Figure 3. The gradient flow coupling on Lja = 14 — 28 volumes. The left panel corresponds 
to ,5 = 1.5, K = 0.14150, the right to /3 = 3.5, k = 0.13380. The strong break between L/a = 24 
and 28 on the left panel indicates that the 28^ volume at these parameter values are in the mass 
deformed regime. 


the t-shift To near its optimal value, tq = 0.15 for the left, tq = 0.1 for the right panel (see 
sec. 4), but setting tq = 0.0 does not change the qualitative results. While on the right 
panel evolves smoothly with the volume, on the left panel we observe a strong break 
between L/a = 24 and 28. This signals that the L/a = 28 configurations have moved into 
the mass-deformed regime and should not be used in the step scaling study. We have seen 
similar behavior at other (5 values. For example at /3 = 2.0, k = 0.13984 we observed that 
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Figure 4. Volume interpolation of in a/L using fourth order polynomials. 

both the L/a = 28 and 24 volumes are in the mass deformed regime. Even at our final 
Kcr = 0.13987 the 24“^ configurations could be mass deformed. Not surprisingly, smaller 
volumes can tolerate larger mistuning of Kcr; and it is also easier to identify an acceptable 
Her value at weaker gauge coupling. Overall, while our k values are not perfectly tuned, 
we can use them on volumes L/a < 24 and /3 > 1.5. Finding on larger volumes and 
stronger couplings would require significant resources and would probably be better done 
using Schrodinger functional boundary conditions. This is beyond our present approach 
and available resources. 

Finally, we mention that we have not found spontaneous chiral symmetry breaking in 
any of our simulations. Even the 28^^, (3 = 1.5, k = 0.14150 configurations that we discussed 
above are deconfined; the gradient flow coupling shown in figure 3 does not show the linear 
increase with the gradient flow time as expected in confining systems. Approaching and 
crossing the Kcr line did not cause numerical difficulties even on the largest volume and 
strongest gauge coupling. We have not observed any first order bulk phase transition 
at strong coupling, neither an Aoki phase nor Sharpe-Singleton type first order phase 
transition [39, 40]. 

4 Analysis of the step scaling function 

For the step scaling function analysis we generated configurations on volumes 8“^, 10^, 
12^, 14^, 16^, 18^, 20^, 24^ and 28^. We found that the 8^ volumes are too small and 
introduce large lattice artifacts, while the 28^ configurations would require more precise 
tuning for (5 < 2.5. We do not use those volumes in our final analysis. We considered 
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Figure 5. The step scaling function at finite L, using volume pairs 10 —>■ 15, 12 —>■ 18, 14 —>■ 21 
and 16 —>■ 24 with c = 0.35, tq = 0.0 and clover operator discretization. The different symbols 
correspond to different (3 values, /3 = 1.5, 2.0, 2.5, .... 6.5 from right to left. The dashed lines are 
fourth order interpolating polynomials in g^. 


11 gauge coupling values, spaced 0.5 apart in the j3 G [1.5, 6.5] interval. We generated 
5,000 — 10,000 thermalized unit-length molecular dynamics trajectories (MDTU) on the 
smaller and 3,000 — 5,000 on the larger volumes, evaluating the gradient flow after every 
10th MDTU. 

In the L/a G [10,24] range we can form 4 volume-pairs with scale change s = 3/2: 
10 —)• 15, 12 —>■ 18, 14 —)• 21, and 16 —)• 24. We do not have L/a = 15 and 21, instead 
we interpolate among the 7 volumes L/a G [10,24]. We considered several interpolating 
functions, polynomials of degree 3 or 4 in either log(L/a) or a/L. Figure 4 shows such 
an interpolation at c = 0.35 and t-shift tq = 0.0, using fourth order polynomials in a/L. 
Our goal is to predict a smooth interpolating form which reduces fluctuations among the 
independent volumes and provide the values of g/. on the missing L/a = 15 and 21. The 
predicted values are independent of the interpolation details, at least within our statistical 
errors. We use these interpolated values to determine the step scaling function /3^i2{g‘l,L) 
of eq. (2.4) at each bare coupling /? and volume pair L —)■ 1.5L. Since the g‘/.{(3,L) values 
are independent, we use standard error propagation to estimate errors. 

Next we perform an other interpolation, this time of I3z/2{9c^ ™ 5c fixed L. We 

consider polynomials of degree 3 or 4 and find that the order of the interpolating polynomial 
does not change the predicted l3z/2{g‘l^L) within statistical errors. An example of such an 
interpolation is shown in figure 5 at c = 0.35 and tq = 0.0, using fourth order polynomials 
for the volume pairs 10 —?■ 15, 12 —?• 18, 14 —?• 21, and 16 —?• 24. Since in the previous 
step we interpolated in L at fixed bare coupling and now we interpolate in g/ at fixed 
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Figure 6. Top panels: {a/L)'^ 0 extrapolation at =4.0, c = 0.35. The left panel shows 

data obtained using the clover while the right panel shows data obtained using the plaquette 
discretization of the energy density operator. Both panels show results at t-shift tq = —0.1, 0.0 and 
0.15. Bottom panels: same as top but for gl =4.8. 


L, there is no correlation and standard error analysis applies. The different symbols in 
figure 5 correspond to different (3 values, (3 = 1.5, 2.0, 2.5, .... 6.5 from right to left. The 
predictions at different /3 are fully independent, yet they show a fairly smooth dependence 
on . The only clear exception is the /? = 2.0 data set, denoted by dashed squares. On 
every volume (3 = 2.0 is off and the deviation increases with the volume. This is a sign of 
mistuning of the hopping parameter k. Any mistiming increases the step scaling function in 
the continuum limit since the mass is a relevant operator. It appears that the step scaling 
function is a better indicator of a non-zero bare mass than the fermion mass discussed in 
section 3. Since in our interpolation we include even clearly mistuned {{3, k) pairs, our 
resulting step scaling function should be considered as an upper limit of the massless case. 

The last step in determining the continuum limit step scaling function is an extrapo¬ 
lation to the (a/L) —)• 0 limit. At fixed g'^ we use the finite volume step scaling functions 
determined in the previous interpolating step. Figure 6 shows this continuum extrapola¬ 
tion at =4.0 and 4.8 with c = 0.35 Wilson fermions can have 0{a) corrections but 
nHYP smearing and the clover term with csw = 1-0 significantly reduce those. We found 

^ We have chosen =4.8 as this is the largest where we can predict /? 3 / 2 (<?ciT) with tq = 0.15 on 
all four volume pairs. As we discuss in section 5.1 the range of g^ decreases with increasing t-shift. 
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Figure 7. The continuum extrapolated c = 0.35 step scaling function using clover discretization, 
To = 0.0, and volume pairs 12 —18, 14 —>■ 21 and 16 —?► 24. The dashed black and magenta lines 
are the 2-loop and 4-loop MS perturbative predictions obtained by integrating the RG /? function 
eq. (4.1). The dotted line is the 2-loop step scaling function as defined in ref. [29]. 


that our values for /3^/2{9c->^) consistent with {a/L)‘^ dependence. The plots show the 
continuum extrapolation at several values of the t-shift parameter, tq = —0.1, 0.0 and 0.15, 
as well as the dependence on the energy density operator (clover versus plaquette). We 
will compare these in section 5. 

Using the clover operator, tq = 0.0, and volume pairs 12 —)> 18, 14 —)• 21 and 16 —)• 24 in 
the range l.S < < 5.2 we obtain the continuum extrapolated step scaling function shown 

in figure 7. For reference we show the universal 2-loop and the MS 4-loop perturbative 
predictions obtained by integrating the continuum RG /3{g'^) function [13] 


r9c+0s(gc) 


dg'^ 


= In(s^), s = 3/2. 


(4.1) 


We also show in figure 7 the 2-loop value obtained by the series given in ref. [29]. The 
more than 15% difference between the two 2-loop discrete /3 functions at g^ = 5.0 show 
that these perturbative predictions should be used as guidance only. From now on we will 
show the perturbative predictions given in eq. (4.1). We will use the results in figure 7 as 
a reference in the comparisons below. 


5 Investigating cnt-off effects 

In this section we investigate the dependence of the continuum limit extrapolated step 
scaling function on the operator used to discretize the energy density, on the gradient flow 
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Figure 8. The continuum limit extrapolated step scaling function Pz/ 2 { 9 t) with tq = —0.1, 0.0 
and 0.15 t-shift, using the clover operator. The top row left and right panels correspond to c = 0.25 
and 0.3, the bottom row left and right panels to c = 0.35 and 0.4. Significant cut-off effects remain 
even after the (a/L)^ 0 extrapolation with c = 0.25, but they are less noticeable with c = 0.3 

and disappear within statistical errors when c > 0.35. The dashed black and magenta lines in all 
panels are the 2-loop and 4-loop MS perturbative predictions. 

defined at different tq t-shift values, and on the volumes used in the (a/L) —)■ 0 extrap¬ 
olation. Any dependence on these quantities indicate that the (a/L) —)■ 0 extrapolation 
did not remove the cut-off corrections and point to higher order 0{a^) or non-perturbative 
effects. Those can be controlled by increasing the parameter c = v^/L, or by using larger 
volumes in the continuum extrapolation. 

5.1 Dependence on the t-shift and the gradient flow c parameter 

Changing tq in eq. (2.2) changes the flow transformation. Deviations in the continuum 
extrapolated values at different tq indicate cut-off effects that are not removed by the 
{a(L)‘^ —)• 0 extrapolation. In section 4 we have discussed the left panels of figure 6 that 
show the continuum extrapolations at g'^ =4.0 and 4.8 using the clover operator with 
Tq = —0.1, 0.0 and 0.15 and c = 0.35. For both g'^ the different t-shifts predict the same 
continuum limit. This observation holds at other couplings as well. The bottom left panel 
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Figure 9. The {a/L) 0 extrapolated step scaling function /33/2(5c) with tq = 0.0, c = 0.35. 

The left panel compares the clover and plaquette discretization, while the right panel compares 
(a/L)^ —>■ 0 extrapolations using all 4 or only the 3 largest volume pairs. The dashed black and 
magenta lines are the 2-loop and 4-loop MS perturbative predictions. 


of figure 8 shows the continuum limit extrapolated step scaling functions with tq = —0.1, 
0.0 and 0.15 and c = 0.35. The agreement between the three gradient flows indicates that 
with c = 0.35 the (a/L)^ —)• 0 extrapolation removed most cut-off effects. It is interesting 
to note that by selecting smaller tq we can reach somewhat larger g'^ . It is possible to push 
this further, selecting tq = —0.25 or even —0.5, and still find consistent results. However 
the slope of the (a/L)^ —?• 0 extrapolation increases with decreasing tq and we do not feel 
comfortable extending the g'^ range this way. 

The other panels of figure 8 show the {a/L)'^ —>■ 0 extrapolated step scaling functions 
with c = 0.25, 0.3 and 0.4 with the same tq = —0.1, 0.0 and 0.15 values. The statistical 
errors are significantly smaller at c = 0.25 on the top left panel. While the dependence of 
/^3/2(5c 5 on (a/L)'^ is still reasonably linear, the iajL'f' —)• 0 extrapolated values depend 
very strongly on the t-shift parameter. Clearly the [ajL)^ —)• 0 extrapolation does not 
remove the cut-off effects at c = 0.25. Had we chosen to work with tq = 0.0 only, we would 
not have observed this inconsistency and could have predicted an incorrect step scaling 
function. 

The step scaling function at c = 0.3, on the top right panel of figure 8, shows much bet¬ 
ter consistency between different t-shift values than c = 0.25. For c = 0.35 and 0.4, shown 
on the bottom two panels, we do not observe any cut-off dependence within our statistical 
uncertainty. Different c values define different RG schemes and their corresponding step 
scaling functions do not have to agree. Nevertheless we see in figure 8 that the continuum 
limit extrapolated step scaling functions with c > 0.3 are very close. 

All panels of figure 8 used the clover operator. The plaquette operator gives very 
similar results. 
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5.2 Dependence on the energy density operator 

The energy density E{t) = —|ReTr in eq. (2.1) can be discretized in many 

different ways. The two most frequently used discretizations are based on the plaquette 
and clover operators [26]. The different discretizations have different cut-off effects and 
comparing the two is an other way to verify that the (a/L)^ —)• 0 extrapolation predicts 
the continuum step scaling function [41]. 

The right panels of figure 6 show the continuum extrapolations at =4.0 and 4.8 
using the plaquette operator and the corresponding C{L,c) values in eq. (2.3). Comparing 
the right and left panels we see that while the cut-off effects at given t-shift are different, 
the continuum extrapolated values are identical with the two operators. This observation 
holds at other g'^ values as well. On the left panel of figure 9 we overlay the continuum 
extrapolated step scaling functions predicted by the plaquette and clover operators with 
c = 0.35, To = 0.0. The agreement between the two operators indicates that with c = 0.35 
the (a/L)^ —)• 0 extrapolation removes most cut-off effects both with clover and plaquette 
discretization. 

Repeating the same comparison at c = 0.25 shows deviation between the two operators, 
but the effect is much less striking than comparing different t-shifts as in figure 8. 

5.3 Dependence on the lattice volume 

Our volumes, while not small, are smaller than ref. [20]. It is important to verify that they 
are large enough to predict the correct continuum limit. The various consistency checks 
discussed above show that with c > 0.35 this is the case. 

In every analysis until now we extrapolated to the continuum limit using volumes 
L/a > 12. As one can infer from figure 6 the L/a = 10 —)• 15 prediction for h/2{9c) 
is fairly consistent with that. (The L/a = 8 ^ 12 predictions, on the other hand, are 
significantly off. We do not show those in figure 6 as we do not use the 8'^ volumes 
anywhere in the analysis.) Including the volume pair L/a = 10 —)• 15 in the {a/L)‘^ —)• 0 
extrapolation indeed has negligible effect on the step scaling function as the right panel of 
figure 9 shows. The agreement is well within statistical errors. 

Looking at figure 5 it is clear that g/. depends only weakly on L at fixed bare coupling 
(5. This is due to the very slow running of the coupling and is independent of the action or 
RG scheme used. The (a/L) —)• 0 limit at larger g/. values on our volumes covers a relatively 
small region in bare couplings, though the interpolation of j3‘^i2{g‘l-,L) in g/ connects the 
weak and strong coupling regions. In any case we should consider the question if the usual 
(a/L)^ 0 extrapolation is acceptable in a slowly running, near-conformal system, or it 

would make more sense to allow an arbitrary power-like dependence in L, as was advocated 
in refs. [16, 33]. While at or near an IRFP one might indeed find a different L dependence, 
we do not see any indication of that in our study. Different t-shift values and operators 
predict the same continuum step scaling function, supporting the validity of the (a/L)^ —)■ 0 
extrapolation. Even in the study of the 12-flavor system an (a/L)^ —)• 0 extrapolation of 
the IRFP is satisfactory [21]. Nevertheless there is a way to mitigate this issue. As figure 6 
shows it is possible to change the (a/L) dependence of /33/2(5c) by changing the tq t-shift 
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value. With the clover operator at c = 0.35 there is practically no dependence on L with 
To = 0.2 in the > 3.5 region, so the actual functional form of the L dependence becomes 
irrelevant. With the plaquette operator this optimal value is also close to tq = 0.2. In 
both cases we find the same continuum extrapolated step scaling function as with other tq 
values. 

6 Conclusion and discussion 

We investigated the step scaling (discrete RG j3) function of the 2-flavor SU(3) sextet 
model, one of the leading candidates of BSM models that describe the Higgs boson as a 
composite particle. Using nHYP smeared clover-Wilson fermions we were able to cover 
the range < 5.5 where we find that the step scaling function follows closely the 4-loop 
(non-universal) MS prediction that predicts an IRFP at g^ ~ 5.90. Lattice calculations can 
suffer from cut-off effects. We made special efforts to investigate them by using different 
discretizations and gradient flow transformations. We found that most cut-off effects are 
removed if the gradient flow parameter c > 0.35, leading to consistent results. 

An independent study of the step scaling function of this model was reported re¬ 
cently [20]. The results, using staggered fermions and the same gradient flow method, 
agree with perturbation theory at small couplings and show a downward deviation relative 
to the 2-loop /3-function at stronger couplings. However this downward deviation is much 
weaker (in fact consistent with the 2-loop prediction of eq. (4.1)) than what we observe 
with Wilson fermions. Since we use the same renormalized coupling to calculate the step 
scaling function as ref. [20] , predictions from different lattice actions are expected to agree. 
However, the results show an over 3a tension at strong couplings between the two works. 

At present we do not understand the source of this disagreement between the staggered 
and Wilson fermion results. This problem should be investigated in the future on larger 
volumes, perhaps using different lattice formulations or even different lattice methods. 

In closing we would like to pose a question regarding universality that might shed light 
to this puzzle. The concept of universality is the cornerstone of critical phenomena and 
lattice QCD calculations, yet it has not been investigated in connection with 4 dimen¬ 
sional many-fermion conformal systems. Universality requires that the lattice actions are 
local, have the same field content, and the same symmetries. If two lattice actions satisfy 
these conditions and are in the vicinity (basin of attraction) of the same fixed point, their 
critical properties (continuum limit) are identical. The two widely used lattice fermion 
formulations, Wilson and staggered fermions, have very different symmetries. The former 
completely breaks chiral symmetry while the latter has a remnant U(l) chiral symmetry. 
Nevertheless in QCD-like systems, around the perturbative Gaussian fixed point, one can 
prove that they differ only in irrelevant operators and describe the same continuum quan¬ 
tum field theory. In conformal systems the fate of the two fermion formulations is not 
obvious in the basin of attraction of the infrared fixed point, as irrelevant operators can 
become relevant and vice versa at different fixed points. 

Similar questions have been studied extensively in 3 dimensional spin systems. A 5- 
loop e-expansion of the 0(ni)©0(n2) model shows that for N = ni+n 2 > 3 the 0{N) fixed 
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point is unstable, the critical properties of the system are controlled either by a biconal 
or a tetra-critical decoupled fixed point [42]. 3 dimensional spin models are very different 
from 4 dimensional gauge-fermion systems, but the richness of the phase structure of the 
3 dimensional spin systems suggests that in 4 dimension we might find unexpected results 
as well. If it turned out that Wilson and staggered fermions have different IRFPs in a 
conformal system, it would not be surprising to find different step scaling function close to 
the IRFP. 

We cannot conclude that the 2-flavor sextet system is conformal, nor that it is chirally 
broken. However our step scaling function approaches zero very closely. If the system is 
not conformal it has to be slowly walking and its IR dynamics could be influenced by a 
nearby IRFP. The question of universality between staggered and Wilson fermions arises 
again. 

Considering the possible phenomenological importance of this system, it is imperative 
to understand and resolve these issues. 
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